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COMPUTED TOMOGRAPHY DEVICE AND METHOD WITH THREE-DIMENSIONAL BACKPROJECTION 



The invention relates to a conaputed tomography method in which an 
examination zone is irradiated along a circular trajectory by a fan-shaped radiation beam and 
in which a detector imit detects primary radiation and radiation which is coherently scattered 
in the examination zone. The invention also relates to a computer tomograph for carrying out 
5 such a method as well as to a computer program for controlling the computer tomograph. 



In known methods of the kind set forth radiation which is coherently scattered 
in the examination zone is measured by a detector unit; the variation in space of the scatter 

10 intensity in the examination zone is then reconstructed from said measuring values. The 
reconstruction is usually carried out by means of iterative methods based on algebraic 
reconstruction techniques (ART) in which a reconstruction of the entire irradiated 
examination zone is necessary. A reconstruction of a part of the examination zone is not 
possible. Moreover, the iterative methods require a large amount of calculation work and 

15 hence lead to comparatively long reconstruction times. Alternative methods involve two- 
dimensional reconstruction techniques based on back projection; in these techniques such 
large distances exist or are assumed between the examination zone and the detector unit that 
an unambiguous relationship exists between the point of incidence of a scattered ray on the 
detector unit and the scatter angle. However, this geometry cannot be realized in medical 

20 applications, so that the two-dimensional back projection produces inaccurate results only. 



Therefore, it is an object of the present invention to provide a method, a 
computer tomograph as well as a computer program which enable a faster reconstruction of 
25 the scatter intensity distribution in the examination zone while taking into accoimt the actual 
geometry in medical applications. 

In respect of the method this object is achieved in accordance with the 
invention by means of a computed tomography method which comprises the steps of 
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- generating, using a radiation source and a diaphragm arrangement which is arranged 
between the examination zone and the radiation source, a fan beam which traverses an 
examination zone or an object present therein, 

generating relative motions, comprising a rotation about an axis of rotation, between 
5 the radiation source on the one side and the examination zone or the object on the 

other side, 

- acquiring measuring values which are dependent on the intensity of the radiation by 
means of a detector imit which detects, during the relative motions, the primary 
radiation from the fan beam and radiation which is coherently scattered in the 

1 0 examination zone or on the object, 

- reconstructing a CT image of the examination zone from the measuring values, during 
which reconstruction a back projection is carried out in a volxmie which is defined by 
two linearly independent vectors of the rotational plane and a wave vector transfer. 

The scatter intensity is dependent not only on the material but also on the 

15 scatter angle and the wavelength of the radiation. A back projection in a volume defined by 
two linearly independent vectors of the rotational plane and a wave vector transfer, however, 
has the advantage that the scatter intensity thus parameterized is dependent exclusively on the 
scatter material, because tiie wave vector transfer is known to be proportional to the product 
of the inverse wavelength and the sine of half the scatter angle. In this context the scatter 

20 angle is the angle which is enclosed by the path of the scattered ray relative to the path that 
the ray would have followed in the absence of the scattering process. 

The scattered rays have a curved shape in said volume. Taking the curved 
shape of the scattered rays into account for the back projection in conformity with claim 2 
enhances the quality of the reconstructed scatter intensity distribution, 

25 According to the preferred reconstruction method disclosed in claim 3 the 

measuring values are weighted prior to the back projection. Such weighting takes into 
account the fact that the effectiveness of the detector surface decreases as the scatter angle 
increases and also that the radiation density decreases as the distance between the point of 
incidence of the scattered ray on the detector unit and the scatter center increases. The 

30 effective detector surface is the detector surface projected in a plane perpendicular to the path 
of the detected ray. The weighting in claim 4 takes into account the fact that the radiation 
density decreases as tiie distance between the radiation source and the scatter center 
increases. The quality of the reconstmcted scatter intensity distribution is enhanced when 
such effects are taken into accoimt. 
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Claim 5 describes a preferred reconstruction method which requkes less 
calcxilation work in comparison with other methods and offers a favorable image quality 
nevertheless. 

A computer tomograph for carrying out the method in accordance with the 
5 invention is disclosed in claim 6. 

Claim 7 defines a computer program for controlling a computer tomograph as 
disclosed in claim 6. 



10 The invention will be described in detail hereinafter on the basis of an 

embodiment and with reference to a drawing; therein: 

Fig. 1 shows a computer tomograph which is suitable for carrying out the 
method in accordance with the invention. 

Fig. 2 is a diagrammatic cross-sectional view of the computer tomograph of 
15 Fig. 1 in the direction of the axis of rotation, 

Fig. 3 shows a flow chart of a version of the method in accordance with the 

invention. 

Fig, 4 is a plan view of a column of detector elements, 

Fig. 5 is a diagranamatic representation of the arrangement of virtual radiation 

20 source. 

Fig. 6 shows a cross-section through the rays of the virtual radiation sources. 
Fig. 7 is a diagrammatic sectional view of the irradiated examination zone and 
the detector in the direction of the plane of rotation, and 

Fig. 8 illustrates the dependence of the amount of the wave vector transfer on 
25 the distance between a scatter center in the exannination zone and the foot of the detector 
(point of incidence of the primary beam on the detector unit). 



The computer tomograph shown in Fig. 1 comprises a gantry 1 which is 
30 capable of rotation about an axis of rotation 14. To this end, the gantry 1 is driven by a motor 
2 at a preferably constant, but adjustable angular speed. A radiation source S, for example, an 
X-ray source, is mounted on the gantry 1. A diaphragm arrangement 31 defines a radiation 
fan beam 41 which is used for the examination and is denoted by solid lines in Fig. 1 . The fan 
beam 41 extends perpendicularly to the axis of rotation 14 and has small dimensions in the 
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direction thereof, for example, 1 mm. Between the diaphragm arrangement 31 and the 
radiation source S there may be arranged a second diaphragm device 32 which forms a 
conical radiation beam 42 jfrom the radiation generated by the X-ray source S. The conical 
radiation beam 42, which would arise in the absence of the diaphragm arrangement 31, is 
5 denoted by dashed lines. 

The fan beam 41 traverses a cylindrical examination zone 13 in which an 
object, for example, a patient on a patient table (both not shown) or also a technical object 
may be situated. After having traversed the examination zone 13, the fan beam 41 is incident 
on a detector unit 16 which is attached to the gantry 1 and comprises a measuring surface 17 

1 0 with a plurality of detector elements arranged in the form of a matrix. The detector elements 
are arranged in rows and columns. The detector columns extend parallel to the axis of 
rotation 14. The detector rows are situated in planes which extend perpendicularly to the axis 
of rotation, that is, preferably on an arc of a circle around the radiation source S. However, 
they may also have a different shape, for example, the shape of an arc of a circle around the 

15 axis of rotation, or be rectilinear. Generally speaking, the detector rows comprise 

substantially more detector elements (for example, 1000) than the detector columns (for 
example, 16). 

The fan beam 41, the examination zone 13 and the detector unit 16 are adapted 

to one another. In a plane perpendicular to the axis of rotation 14 the dimensions of the fan 
20 beam 41 are chosen to be such that the examination zone 13 is fully irradiated and the length 

of the rows of the detector unit 16 is such that the fan beam 41 can be detected completely. 

The fan beam is incident on the central detector row (rows). 

When a technical object is concerned instead of a patient, the object can be 

rotated during an examination while the radiation source S and the detector unit 16 remain 
25 stationary. The examination zone 13, or the object or the patient table, can be displaced 

parallel to the axis of rotation 14 by means of a motor 5. As an equivalent, however, the 

gantry could also be displaced in this direction^ 

When the motors 2 and 5 operate simultaneously, the radiation source S and 

the detector unit 16 describe a helical trajectory relative to the examination zone 13. 
30 However, when the motor 5 for the displacement in the direction of the axis of rotation 14 is 

stationary and the motor 2 rotates the gantry, a circular trajectory is obtained for the radiation 

source S and the detector 16 relative to the examination zone 13. Hereinafter only the circular 

trajectory will be considered. 
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Fig. 2 shows that between the examination zone 13 and the detector xinit 16 
there is arranged a collimator arrangement 6 which comprises a plurality of blades 60. The 
blades 60 consist of a material which strongly absorbs the X-rays and are situated in planes 
which extend parallel to the axis of rotation 14 and intersect one another at the focus of the 
radiation soxirce S. Their spacing may amount to, for example, 1 cm and each blade 60 may 
have a dimension of, for example, 20 cm in the plane of drawing. The collimator arrangement 
6 thus subdivides the fan beam 41 into a number of neighboring segments, so that essentially 
a column of detector elements is struck by primary or scattered radiation from one segment 
only. 

The measuring values acquired by the detector unit 16 are applied to an image 
processing computer 10 which is connected to the detector unit 16, for example, via a 
wireless data transmission. The image processing computer 10 reconstructs the scatter 
intensity distribution in the examination zone 13 and reproduces it, for example, on a monitor 
11. The two motors 2 and 5, the image processing computer 10, the radiation source S and 
the transfer of the measuring values from the detector unit 16 to the image processing 
computer 10 are controlled by a control unit 7. 

In other embodiments the acquired measuring values can first be applied to 
one or more reconstruction computers for reconstruction, which reconstruction computers 
apply the reconstructed data to the image processing computer, for example, via an optical 
fiber cable. 

Fig. 3 shows a flow chart illustrating the execution of a version of a measuring 
and reconstruction method that can be carried out by means of the computer tomograph 
shown in Fig. 1. 

After the initialization in the step 101, the gantry rotates at a constant angular 
speed. In the step 103 the radiation of the radiation source S is switched on so that the 
detector unit 16 can detect tibie prLtnary and the scattered radiation from a plurality of angular 
positions. The detector element or elements at the center of each detector colunm detect 
essentially the primary radiation, whereas the scattered radiation (secondary radiation) is 
detected by the detector elements which are situated further nearer to the ends of a column. 

The foregoing is diagrammatically indicated in Fig. 4 which is a plan view of a 
column of detector elements. The detector elements 161 which detect the scattered radiation 
are denoted by simple shading whereas the detector element 160 which is situated at the 
center and detects the primary radiation is marked by a cross. In other radiation sources, 
notably in radiation sources having a larger focus, the primary radiation could also be 
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detected by more than one detector element To both sides of the central detector element 
there are situated detector elements which are struck by scattered radiation, but also by a 
(reduced) amount of primary radiation (white elements neighboring 160) because of the finite 
dimensions of the focus of the radiation source. In this embodiment scattered radiation is 
5 considered to be only the radiation which is measured by the detector elements denoted by 
shading in the drawing. 

The scatter intensity is dependent inter alia on the energy of the scattered X- 
ray quantum. Therefore, the energy of the scattered X-ray quantum must either be measured, 
implying that the detector elements should be capable of measuring in an energy-resolved 
10 fashion, or use must be made of X-rays with quantum energies fi:om an as small as possible 
range (monochromatic X-rays in the ideal case). There are various possibilities for 
niiiiimizing the energy difference of the X-ray quanta relative to their energy: 

- The use of suitable filter materials, for example, copper in the primary beam. As a 
result, the soft X-rays generated by an X-ray source, that is. X-rays having a low 

1 5 quantum energy, are suppressed to a high degree. 

- Additionally, the voltage of an X-ray tube can be optimized with respect to the 
selected filter. 

Finally, the so-called "balanced filter" technique can be applied. According to this 
technique the data is acquired twice and each time filters which have a slightly 
20 different atomic number and whose K edge is used for filtering are then situated in the 

beam path. Subsequently, the difference signal firom the two measurements is 
extreicted. 

In tiie step 105 the measuring values of the scattered radiation are normalized. 
The measuring values of each radiation source position of the scattered radiation are then 

25 divided by the measuring values of those primary rays which have caused the scattered rays. 

In the step 107 the measuring values are subjected to a rebinning operation. 
Each measuring value is Ihen associated with a line firom the detector element in which the 
measuring value has been detected to the radiation source position. Thus, it is assumed that 
rays of the fictitious, conical radiation beam 42 have caused the measuring values without the 

30 rays having been scattered. Due to the rebinning operation the measuring values are resorted 
£is if they had been measured by means of a different radiation source (a circular radiation 
source capable of emitting mutually parallel fan beams) and a different detector (a flat, 
rectangular virtual detector). 
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This will be described in detail with reference to Fig. 5. Therein, the reference 
numeral 17 denotes the circular trajectory along which the radiation soiurce irradiates the 
examination zone. The reference numeral 413 denotes a fan-shaped radiation beam which 
emanates from the radiation source position So and whose rays propagate in a plane 

5 containing the axis of rotation 14. The conical radiation beam emitted by the radiation source 
in the position So may be assumed to consist of a pluredity of flat fan beams which are 
situated in planes which are parallel to the axis of rotation 14 which intersect in the radiation 
source position. Fig. 5 shows only a single one of these fan beams, that is, the fan beam 413. 

Moreover, Fig. 5 also shows further fan beams 41 1, 412 and 414, 415 which 

10 extend parallel to the fan beam 413 and are situated in planes which are parallel to one 

another and to the axis of rotation 14. The associated radiation source positions S-2, S-i and 
Si, S2 are occupied by the radiation source S before and after reaching the radiation source 
position So, respectively. All rays in the fan beams 41 1 to 415 have the same projection 
angle. The projection angle is the angle enclosed by the planes of Ihe fan beams with respect 

15 to a reference plane which is parallel to the axis of rotation 14. 

The fan beams 41 1 to 415 define a radiation beam 410 having a tent-like 
shape. Figs. 5 and 6 show the sectional area 420 which is obtained when the radiation beam 
410 is intersected by a plane which contains the axis of rotation 14 and extends 
perpendicularly to the planes of the fan beams 41 1 to 415. The two edges of the sectional 

20 area 420 which intersect the axis of rotation are curved. This curvature is due to the fact that 
the radiation source position at the center (for example. So) are situated further from the 
sectional plane than those at the edge (for example, S2 or S.2) and that the fan beams all have 
the same angle of aperture. For each group of fan beams, therefore, a rectangular virtual 
detector 170 is defined in the plane sectional area 420, the edges 171 and 172 of said virtual 

25 detector being given by the dimensions of the outer fan beams 411 and 41 5 in the plane 
sectional area. 

Fig. 6 also shows (marked by round dots) the puncture points of some rays, 
contained in the fan beams 41 1 ... 415, through said virtual detector. Finally, crosses mark the 
supporting points of a regular cartesian grid. The puncture points and the supporting points 
30 usually do not coincide. From the measuring values for the puncture points, therefore, the 
measviring values at the eqmdistant supporting points within the virtual detector 170 have to 
be determined. Such rebinning is described in detail in DE 198 45 133 Al whereto reference 
is explicitly made herein. 
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In the step 1 09 the measuring values of the scattered rays as produced by the 
rebinning are subjected to one-dimensional filtering with a factor which increases ramp-like 
as a function of the spatial frequency. To this end, each time successive values are taken into 
account in a direction parallel to the plane of rotation, that is, along a row of the virtual 
5 detector. This filtering is carried out along each row of the virtual detector for all projection 
angles. 

In other embodiments the rebinning could also be omitted. It is knoAvn that in 
that case the filtering must be modified because the detector imit is curved, for example in the 
form of an arc aroimd the radiation source or around the axis of rotation. 

1 0 The filtered measuring values are subsequently xised for the reconstruction of 

the scatter intensity distribution in the examination zone by back projection. 

The back projection is carried out in a volume which is defined by the vectors 
X , y and q , the unity vectors x and y being situated in the rotational plane and being 
oriented horizontally and vertically, respectively, whereas the wave vector transfer g is 

15 oriented parallel to the axis of rotation. In other embodiments two other, linearly independent 
vectors of the plane of rotation can be used instead of the vectors x and y . As has already 
been stated, the amount of the wave vector transfer q is proportional to the product of the 
inverse wavelength X of the scattered X-ray quanta and the sine of one half of the scatter 
angle ©: 

20 q = (l/X)sin(0/2) (1) 

The scatter angle © can be determined on the basis of the arrangement of the 
examination zone 15 irradiated by the fan beam 41 and the detector imit 16 as shown in Fig. 
7. As opposed to known two-dimensional methods based on back projection, this Figure 
takes into accoimt the real geometry in medical applications. For one or each detector 
25 element Dj this gives rise to different scatter angles © which are calculated in conformity 
with the following equation: 

0 = arctan(a/d) (2) 
Therein, d denotes the distance of a scatter center Si and a denotes the distance between the 
detector element D\ and the foot 12 of the detector. 
30 The detector Di detects rays which have been scattered at the angles 

©1 < 0 < ©2 in the examination zone 15 irradiated by the fan beam 41 . 

From the above two equations the following equation is obtained for small 

angles © 
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q - a/(2dX) (3) 
The variation of the value of the wave vector transfer q in dependence on the 
distance d between a scatter center and the foot of the detector thus has the shape of a 
hyperbola and hence is not linear (shown in Fig. 8). It follows therefrom that the originally 
5 linear path of the scattered rays in the (x,y,q) space is curved. The back projection thus takes 
place along rays curved in the form of a hyperbola. 

In the step 111a voxel V(x,y,q) is determined within a selectable (x,y) zone 
(field of view or FOV) and within a range of values of the wave vector transfer q resulting 
from the geometry of the computer tomograph. 

10 In the step 113 the filtered values are multiplied by a weighting factor which 

corresponds to the reciprocal value of the cosine of the scatter angle. The decreasing effective 
detector surface as a ftinction of an increasing scatter angle is thus taken into account. When 
said angle is small, the cosine of the angle is practically always 1, so that this weighting can 
be dispensed with. Moreover, the fact that the radiation density decreases as the distance 

15 between the scatter center increases, that is, the voxel V(x,y,q), and the point of incidence of 
the scattered ray on the measuring surface increases is taken into accoimt by multiplying all 
measuring values for each radiation source position by a weighting factor which corresponds 
to the square of the distance between the scatter center, at which the ray associated with the 
measuring value was scattered, and the point of incidence of the scattered ray on the 

20 measuring surface. 

When the rebinning is omitted in other versions, an additional multiplication 
of the filtered measuring values by a weighting factor is necessary; this factor is equal to the 
reciprocal value of the square of the distance between the radiation source position and the 
scatter center at which the detected ray was scattered. 

25 During the back projection in the step 11 5 all curved rays are then taken into 

accoxmt which pass through the voxel V(x,y,q). If not one ray from a radiation source 
position passes exactly through the center of the voxel, the associated value must be 
determined by interpolation of the measuring values of neighboring rays. The measuring 
value that can be associated with the ray passing through the voxel, or the measuring value 

30 obtained by interpolation, is accumulated on the voxel V(x,y,q). After the contributions to the 
relevant voxel have thus been accumulated for all radiation source positions, it is checked in 
the step 1 17 whether all voxels in the (x,y,q) zone to be reconstructed have been dealt with. If 
this is not the case, the flow chart branches to the step 111. Otherwise the scatter intensity 
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distribution has been determined for all voxels in the FOV and the reconstruction method is 
terminated (step 1 19). 

In other versions the back projection in the volume defined by two linearly 
independent vectors of the rotational plane and the wave vector transfer can be carried out as 
5 an approximation along straight rays. 



